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'Zl Abstract 



We develop a complete mathematical theory for the symmetrical solutions of the 

generalized nonlinear Schrodinger equation based on the new concept of angular 

pseudomomentum. We consider the symmetric solitons of a generalized nonlinear 

Schrodinger equation with a nonlinearity depending on the modulus of the field. 

We provide a rigorous proof of a set of mathematical results justifying that these 

solitons can be classified according to the irreducible representations of a discrete 

group. Then we extend this theory to non-stationary solutions and study the re- 

lationship between angular momentum and pseudomomentum. We illustrate these 

theoretical results with numerical examples. Finally, we explore the possibilities of 

J> the generalization of the previous framework to the quantum limit. 
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1 Introduction 



One-dimensional discrete solitons in periodic media have been proposed both 
in the framework of nonlinear optics [Tj and that of Bose-Einstein conden- 
sates (BECs) PIB] . On the other hand, their two-dimensional counterparts 
have attracted intense investigational effort in recent years. In this context, 
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it has been shown that two-dimensional optical discrete solitons can exist in 
optically induced lattices [HE], in nonlinear photonic crystals piT] . and in 
nonlinear photonic crystal fibers [S]. In addition, experimental observations 
of these solitons in one and two dimensions have been reported in optically 
induced lattices pilOIIll] and in modulated Bessel optical lattices [T2IT3] . In 
parallel to these findings, discrete soliton vortices have been proposed both in 
the context of discrete models |T3] and in continuous models of optically in- 
duced lattices p[T3] , Bessel lattices [10] , and photonic crystal fibers [17j . They 
have been also experimentally observed in optically induced lattices p^l9f20] . 
Other discrete solitons, such as dipole and quadropole solitons, have been 
shown to exist by means of numerical simulations in optically induced lat- 
tices [T5|21f22j . in Bessel lattices [23], and in photonic crystal fibers [23], and 
they have been observed in the former system [21j- Bright discrete solitons, 
vortex solitons and dipole solitons have been also proposed in a self-attractive 
(or focusing) BEG loaded in an optical lattice [5)1151125] . Different results have 
been obtained in the self-repulsive (or defocusing) case. Extended solitons 
or BEG arrays have been realized in 2D optical lattices [26|27] . These BEG 
arrays displayed coherent spatial structures in 2D induced by modulational 
instability [28j. This was the mechanism that permitted researchers to pro- 
pose in-phase and vr out-of-phase discrete bright solitons in self-defocusing 
media [29130]. 

All these localized structures have been found as soliton solutions of 2D non- 
linear Schrodinger Equations (NLSE) with a periodic potential and different 
expressions for the nonlinear part (Kerr focusing or defocusing nonlinearity, 
saturable, etc) - we will refer to all such equations as generalized nonlinear 
Schrodinger Equations. In addition, the amplitude of all these solitons has dis- 
crete rotational symmetry. In Ref. [23] it was proven that a soliton of a gener- 
alized NLSE in a discrete rotationally symmetric potential with a nonlinearity 
depending on the modulus of the field must belong to some representation of 
the corresponding discrete rotational group, provided that its modulus is also 
discrete rotational invariant. As proved in Ref. [31], this property implies that 
there must exist a "cut-off" in the allowed values of the vorticity of a discrete 
vortex soliton in discrete rotational symmetry media, as long as the discrete 
vortex presents a single singularity. It was also demonstrated that identical 
results can be obtained if the discrete vortex soliton is reinterpreted as an 
angular Bloch mode characterized by its angular pseudomomentum, which is 
conserved along the direction of propagation ^3]. Another important result 
is the existence of a rule that links angular momentum and angular pseu- 
domomentum in an interface between a rotationally invariant and a discrete 
rotationally invariant medium [32]. This rule permits the study of vorticity 
transformations in this kind of interfaces [3 21IM] and some of these transfor- 
mations have been experimentally observed [35] . 

In this paper we will rigorously prove that the solitons of a generalized NLSE 
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with a nonlinearity depending on the modulus of the field can be classified 
according to the irreducible representations of a discrete group, if that mod- 
ulus is discrete rotationally invariant. We will show that, consequently, a dif- 
ferent quantity can be assigned to discrete solitons, discrete vortex solitons, 
and discrete /c-pole solitons (for example, dipole or quadrupole solitons). We 
will name this quantity angular pseudomomentum to be consistent with the 
nomenclature of Ref. [33]. We will also offer detailed proofs of the conserva- 
tion of angular pseudomomentum and transformation rules between angular 
momentum and pseudomomentum. Finally, we will outline the extension of 
these results to the quantum limit. 

The paper is organized as follows. In the Sec. |2j the classification of discrete 
solitons according to angular pseudomomentum is demonstrated. In Sec. [3] 
examples of this classification are provided. Section |4] is devoted to the conser- 
vation of angular pseudomomentum and to the transformation rules between 
angular momentum and pseudomomentum. In Sec. [5] the extension of the pre- 
vious results to the quantum limit is outlined. Finally, in Sec.[6]the conclusions 
are presented. 



2 Classification of Discrete Rotational Symmetric Solitons 



Let Lo(x) = A - F(x), where x = (x, G M^ A = ^ + |J and \/(x) 



is a function invariant under a symmetry group Q. Let us assume that Q is 
the rotational symmetry group 0{2) (not to be confused with an order of 
magnitude symbol), a discrete rotational symmetry group C„ of order n, or a 
discrete rotation- refiection symmetry group Cnv- 

Let Li{\ip\'^) = Fdipl"^), where ip = ipix.) is a function such that IV'p is invariant 
under a symmetry group Q' and F is a function invariant under the same group 
as I'?/' p. We will assume that Q', the symmetry group of l^pul"^, is a subset or 
coincides with Q, the symmetry group of ^(x). Let L(x, l^/'p) = Lq + Li. 
Under these assumptions, let us consider the eigenvalue equation 



A solution ip of this equation gives rise to a stationary solution 0(x, t) = 
e*^V(x) of the generalized NLSE 



Let ipui'^) be a solution of Eq. ([T]). We call the differential operator gener- 
ated by ipy the function L(x, |'?/'i/p). The eigenfunctions ifj of this operator 




(1) 




(2) 
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with eigenvalue /i satisfy L(x, = nip. The function V'i^(x) is a se//- 

consistent solution since L(x, l'?/';^^)'?/',^ = ^ipu- The function ■i/',^ is a symmetric 
self- consistent solution if ^' C ^. A function -0^ is a non- self- consistent solu- 
tion if L(x, |'0iyp)^/'T; = A^V'r? but L(x, I'j/'^p)?/'^ 7^ /i?/^^. 

Note that the differential operator L(x, l^/'p) is Hermitian since L = L\ and 
therefore all its eigenvalues /i are real, G M. Moreover, L = L* since l^(x) 
and Fd'i/'p) are real functions. 

For a self-consistent solution ip^, if we conjugate the eigenvalue equation ([T]) 
we find that Lq'?/'* + I'?/'* P )'?/'* = since is a real function and {ipul"^ = 

F is a real function, and therefore, Lid^p^l"^) = Li(|'?/'*p). Hence, ipt is 
also a self-consistent solution with eigenvalue /i* of the same operator. Finally, 
fi = fi* due to the hermiticity of L. This proves the following Lemma: 

Lemma 1 // ipi, is a self-consistent solution with eigenvalue fi, ipl is also a 
self- consistent solution with the same eigenvalue. 

Let us prove now the following important Lemma: 

Lemma 2 // ip^, is a self-consistent solution with eigenvalue fi, and iprj is a 
function such that L(x, \ipij\^)ipfj = fia'^'a, then ip^j is a self- consistent solution 
if and only if (i) ipy = kip^ or, (ii) ip^ = kip*, with k E C and = 1. 

If ipa is a self-consistent solution then 

L(x, \ip^\'^)ip„ = fiip„. 

By hypothesis 

Now we substract the equations above to obtain 
F{\iP,\')iP^ - F{\iP^\^)iP^ = 0, Vx. 

For most physical systems F is an injective function (i.e., one to one) and 
therefore \ipu\'^ = iV'o-P, Vx. This fact implies (i) ipi, = kipa or, (ii) ipu = kip* 
with /c G C and \k\ = 1. 

On the other hand, if (i) ipi, = kip„ or, (ii) ip^ = kip^., obviously ipa- is a self- 
consistent solution since, in both cases, L(x, \ip„\'^)ip„ = fiip„. This proves the 
previous Lemma. 

Now we proceed to apply this Lemma to the transformed functions T/ip^,, 
where T/ is an element of the group Q' . Before doing this we observe that 

7;'L(x, l^^n^^(x) = 7;V^V^^(x), 
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or 



7:'L(x, |V^,n7:'-i7-V.(x) = /i.7:V.(x). 

It is clear that L(x, \ipu\'^) is invariant under the symmetry group Q' of {ipiyl"^- 
Firstly, A is invariant under the symmetry group 0{2), and consequently it is 
invariant under any discrete rotational symmetry group. Secondly, the function 
V{x.) is invariant under the symmetry group Q. Since we have assumed that 
Q' G Q it must be also invariant under the symmetry group Q'. And finally, 
Fdip^l"^) is invariant under the symmetry group Q' since {ipul"^ is also invariant 
under this group. So T/L^x, \il)y\'^)Tl~^ = L(x, lipu]"^) and therefore 

L(x, \i;,\')T:M^) = f^VM^)- 

Hence T/ilJiy{x.) is also an eigenfunction with the same eigenvalue. Then, the 
conditions of Lemma |2] are satisfied for -ipa- = T/'tp^{'x). 

We conclude that ipa is a self-consistent solution verifying 

if and only if ipa = %'4'u = kipu or ip„ = %'ipu = kipt with \k\ = 1. 

Proposition 3 Let ip^ix) be a symmetric self- consistent solution of Eq. 
Then, for any element T/ of the group Q' , the function T/ipi, can be only a 
self- consistent solution of Eq. ([7]j if and only if it is proportional to or to 
ip*, up to a complex number of modulus unity. 

We have used that Tj'L{x, \-ip,,\'^)T-~^ = L(x, iV^i/P) in the demonstration of 
the previous proposition. This is equivalent to the condition 

[l,t:] = o. 

It is well known in quantum mechanics that if %' and L commute, we can have 
simultaneous eigenfunctions of both operators. This is equivalent to the state- 
ment that the eigenfunctions of L characterized by the eigenvalue /i belong to 
the irreducible representations of the group ^'[36]- The representations of Cn 
or Cnv are one or two-dimensionaQ In the case of two-dimensional represen- 

^ All the irreducible representations of a discrete rotational group C„ are one- 
dimensional. Since L = L* , there are pairs of irreducible representations whose 
characters are mutually conjugate and that share the same eigenvalue [36]. In quan- 
tum mechanics this is due to the time reversal invariance of the Hamiltonian [36] . 
We will consider these pairs of irreducible representations as a doubly-degenerate 
representation. 
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tations, the rotational elements of these groups are given by two-dimensional 
unitary matrices. In the diagonal basis this matrices are of the form 



where e = e*^. On the other hand, we have shown that the self-consistent 
solution tpj, verifies that T/ipu = kip^, where \k\ = 1. Therefore, ipi, is an eigen- 
function of T/ with eigenvalue k. Comparison with the previous expression 
implies that k = e. Moreover, group theory insures that the other eigenfunc- 
tion of the operator has to be ipl [36j- We can summarize these arguments 
in the following Theorem: 

Theorem 4 Let ipv = V'i'(x) be a symmetric self-consistent solution and let us 
assume that Q' <ZQ . Then ip^ belongs to some of the irreducible representations 
ofQ'. 

But, moreover, we have proved also the following Theorem: 

Theorem 5 Under the conditions of the previous Theorem and for two-dimen- 
sional representations, the only self- consistent solutions are those that form 
the diagonal basis, i.e., the pair ipi, and ipl that fulfills, on the one hand, 
L(x, \ipu\'^)i^u = fJ'i^u and L{-x, \ipl\'^)ipl = fiipl, and, on the other hand, Ti'ipu = 
eip^ and T/ipl = e*i)l. 

Remark 6 Note that if Q' d Q the symmetric self- consistent solutions belong 
to irreducible representations of the subgroup Q' . This is the case of the nodal 
solutions introduced in Ref. [24^ or the dipole solutions introduced in 
Refs. IMMMM- 

Remark 7 Discrete rotation-reflection groups C^v present one- and two-dimen- 
sional representations. Under the p^^ rotational element of the group the 
function transforms as C^ip = kpip where kp E C is a unitary constant. But 
under the reflection elements of the group U^, the function is transformed as 
Huip = kjjip* or as H^'ip = kyip, where k^ E C is a unitary constant. So these 
transformations could involve the function and its complex conjugate. 

Discrete rotational groups, C„, present | -l- 1 irreducible representations for 
even n and + 1 for odd n, if the mutually conjugated representations can 
be considered as a single two-dimensional representation. Two of these repre- 
sentations are one-dimensional, while the rest are two-dimensional. Adding re- 
flections to Cn to form a Cnv group means, on the one hand, that these mutually 
conjugated representations form a single two-dimensional representation. On 
the other hand, the two one-dimensional representations split due to different 
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behaviours of the functions with respect to the reflection axes of the group [36] . 
In both kinds of groups, all the transformation rules can be satisfied by func- 
tions of the form il){r, 9) = e*^'"'u(r, 6), where u{r, 9) = u{r, 0+'^) and (r, 6) are 
polar coordinates. Therefore, u{r, 9) is some periodic function in 9 with period 
^. Moreover m G {0, ±1, ±2, . . . , f } for even or m G {0, ±1, ±2, . . . , ±^} 
for odd n according to Th. |8j Functions with m = belong to the first one- 
dimensional representation of C„ or to some of the one-dimensional represen- 
tations of Cnv Each of these one-dimensional representations have a different 
behaviour under reflections. Functions with m = |, n even, belong to the other 
one-dimensional representation of C„ or to some of the one- dimensional rep- 
resentations of Cnv The rest of the functions with m G {±1, ±2, . . . , ±| =]= 1} 
for even n or {±1, ±2, . . . , ±^^} for odd n belong to two-dimensional repre- 
sentations. This can be summarized in the following Theorem: 

Theorem 8 (Cut-off for m). The modulus of m presents a maximum value, 
i.e., \m\ < I for even n and \m\ < for odd n. 

We will name the variable m angular pseudo-momentum, in accordance with [53] 
We will call solutions with m = 0, fundamental solitons. We will call solutions 
with m = |, for even n, nodal or fc-pole solitons. Finally, we will call the 
solutions with m equal to one of the other allowed values vortex solitons. 



3 Classification of Solutions According to their Angular Pseudo- 
momentum in C^-u and Cg^ 

Let us numerically calculate some solutions in C4t, and Cq^ with all possible an- 
gular pseudomomenta using a relaxation numerical method. We will calculate 
only solutions with no more than one singularity and with trivial behaviour 
under reflections. We leave the study of solutions with more singularities and 
different behaviour under reflections to future work. We will solve Eq. ([T]) with 
F{\ip\'^) = 5f(x)|'?/'p, where (^(x) = K(x), i.e., in the self-focusing case. 

3. 1 Solutions in C^y 

In this case ^(x) = Vq [cos^(x) + cos^(?/)], with Vq = 2. This potential can 
model the propagation of light in an optically induced lattice over a photore- 
fractive crystal or a BEC loaded in an optical lattice in the attractive case. 
In both cases appropriate transformations in the physical variables must be 
done to obtain Eq. (bj). We will consider solutions with invariant under 
the same group as V^(x), i.e., C^y. Therefore, n = 4 and the allowed values of 
m according to Th. |8]are m G {0, ±1, 2}. In Fig. [ijthe amplitude and phase of 
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(b) 




Fig. 1. Symmetric self-consistent solutions when Q and Q' are €4^, and /x = 1. 
The amplitude of the solution with angular pseudomomentum m E {0,-1,2} is 
represented in (a), (c), and (e), respectively. In (b), (d), and (f) the corresponding 
phases of these solutions are represented. 

these solutions are presented. They are described by the character table of the 
corresponding irreducible representation since they are correctly transformed, 
i.e., as e^"^'^, under the action of the rotational elements of the group, i.e., 
rotations of ip = radians with p G Z. The fundamental soliton (m = 0) 
and nodal soliton (m = 2) solutions belong to one-dimensional representa- 
tions. The vortex soliton solution (m = —1) belongs to a two-dimensional 
irreducible representation; its complex conjugate is also a solution, but with 
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m = 1. 



3.2 Solutions in Cq^ 

In this case V^(x) is equal to one in circles of radius R = j, where R = 8000 
and A = 1064, and zero outside of these circles. The centers of the circles are 
located in an hexagonal grid, i.e., in positions r = a[kx + cos{a)]i + aky sin(a)j, 
where i, j define a cartesian coordinate system, fc^, /cy G Z, a = |, and a = 
with A = 23000. This ^(x) is used to model the propagation of light beams 
in photonic crystal fibers as in Refs. Pll7|24|31j . 

We will consider solutions with l^pf invariant under the same group than 
V(x), i.e., Cqv For n = 6 the allowed values of m acording to Th. |8] are 
m e {0, ±1,±2,3}. Fig. [2] shows the amplitude and phase of self-consistent 
solutions with m G {0, —1, —2, 3}. It is easy to see that these solutions fulfill 
the character table of the corresponding irreducible representation using any 
rotation of rad with p G Z. For m = or m = 3 the solutions are 
fundamental and nodal solitons, respectively. That is, they belong to one- 
dimensional representations of the C^y group. For m G { — 1, —2} the vortex 
soliton solutions belong to two-dimensional irreducible representations. Their 
complex conjugates are also solutions but with m G {1,2}, respectively. 



4 Dynamics of Rotationally Symmetric Solutions 

4-1 Conservation of angular pseudomomentum 

Let 0(x, to) be a stationary or non-stationary solution of Eq. Q in certain 
time to- Let its modulus |0(x, to)!^ be invariant under Cn or C„^, and let us call 
its symmetry group Q' . As usual, we assume that this group coincides with 
or is a subgroup of the symmetry group Q of V^(x), where Q can be C(2), C„, 
or Cnv Using similar arguments as the ones presented in the previous section 
it can be shown that such a solution must have angular pseudomomentum m 
and that the value of m must obey Th. |8j Therefore, these solutions are of the 
form 0(x, to) = e*™^u(r, 9, to), with u{r, 6 + e,t) = u{r, 6, t) and e = ^• 

If 0(x, to) is a stationary solution, angular pseudomomentum is trivially con- 
served for t > to. Let us prove that if it is not a stationary solution the angular 
pseudomomentum is also conserved. We will prove this statement by induc- 
tion. Let us denote 0o by 0(x, to). It can be shown, using the same arguments 
as in Sec. pi that the operator L(x, |0o|^) = Lq + Li is invariant under the 
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(a) 



(b) 





-pi/3 
-2pi/3 



(d) 




(f) 




(g) 




(h) 



-pi' 




Fig. 2. Symmetric self-consistent solutions when Q and Q' are Cq^, and /i = —0.06. 
The amplitudes of the solutions with angular pseudomomentum m £ {0, —1, —2, 3} 
are represented in (a), (c), (e), and (g), respectively. The circles in which V{x.) is 
equal to one are represented in light blue. In (b), (d), (f), and (h) the corresponding 
phases of these solutions are represented. 
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rotational elements of the group Q', which must be C„ or Cnv Therefore, 
C^LC^^ = L, where = radians, p E Z. 

Let us discretize the evolution from to = to t using arbitrarily small steps of 
size h with h ^ 0. Let tj = to +j h, j G {0, 1, . . .}, and (pj = 0(x, tj). In a first 
order equation in the evolution variable, the solution (pj+i can be obtained in 
terms of (f)j as = e'^^'^^'^'^^^ ^^4>j- 

On the one hand, C^L(x, \M^)C^^ = L(x, \(j)of) and therefore C^e^-^^'^'I'^ol'^'^C^^ = 
giL(x,|0o| since for h ^ and using the series expansion for the expo- 
nential (l + iL{-x, |0o|^)^) = 1 + iL(x, \(j)Qf)h. Consequently, C^0i = 
(-^giL{x,|<^oP)/^0g = e*-^(^'l'^ol')'^C^0o. Taking into account that C^^o = e'"'^(j)o it 
follows that Ci)(f)i = e^^(^'l<^ol')?«eim,?0^_ Consequently, C^^i = e'"''^(f)i. Then 
is a function with the same angular pseudomomentum as 0o- 

On the other hand, let us suppose that (pj is a function with well defined 
angular pseudomomentum m. Then C^4'j = e^"^^(j)j and C^L(x, |0j|^)C^^ = 
L(x, Using the discrete evolution equation and these two equations in 

the same way as in the previous paragraph, one finds that C^0j+i = e*'"''0j+i 
and therefore (pj+i and (pj share the same angular pseudomomentum. In the 
continous limit h —>■ both results are true. This proves the following Theo- 
rem: 

Theorem 9 Let 0(x, t) be a solution of Eq. ^ and let us assume that Q is 
0(2), Cn or Cnv Let us assume as well that the symmetry group Q' of the 
solution at certain to, 0(x, tg); is Cn or Cnv one? Q' C Q. Then the solution 
presents angular pseudomomentum m and is conserved for t > to. 

A solution with well defined angular momentum must be an eigenfunction 
of TZ(f) = e*'^0, where / is the angular momentum and TZ = e~^^^ae\ If the 
evolution operator L is invariant under 0(2), I is conserved and coincides 
with the z component of the angular momentum: 

3z = -^J dxdy^*{r x V)(^ , (3) 

R2 



where (/> is a normalized eigenfunction of TZ. In a system described by a nonlin- 
ear evolution operator L as defined before, this is no longer true if the function 
0(x, to) is such that |0(x, to)P is invariant under C„ or Cnv, even in the case 
where Q is 0{2). Of course, it is not true if Q is C„ or Cnv 
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4-2 Transformation rule 



Let us study the solutions of Eq. ^ when the solution has a well defined 
angular momentum for t < ti and the symmetry of V(x, t) changes into a 
discrete group C„ or Cnv at ti, i.e.: 

l/(r,^,t)= ^ " (4) 

\/(r,^), I V{r,e) = V{r,e+^-^) ift>ti. 



The function 0(x, t < ti) has well defined angular momentum and ^(x, t = ti) 
also has well defined angular pseudomomentum m, since any discrete rota- 
tional group is a subgroup of 0{2). Therefore, C^0(x, z = zi) = e™'^0(x, t = 
ti) with = p"^ radians, where p G Z and is a rotation element of C„ or 
Cnv This m is conserved for t > ti due to Th. [9j Let us see that the final value 
of m depends on the initial value of I and on the order of symmetry n. Let 
4'i{'^,t < ti) be the function with angular momentum / and (/)m(x, t > ti) the 
function with angular pseudomomentum m. The projection coefficient of (pi on 
is = ^i)0/(x, tOd^x. Also, 0„(r, + ^, ti) = e^'"^0™(r, ^, ti) 

" (f)i{r,6,ti). So if we take the transformation 6 
6* + ^ in the expression of Cmi we obtain c^i = e**^'~™'^^Cm;. Therefore the 
coefficient must be zero for all m except those that fulfill I — m = kn with 
G Z. On the other hand, m must present a cut-off, in accordance with Th. 
[8| This proves the following proposition: 

Proposition 10 Let us suppose that V{'x, t) is of the form given by Eq. 
and that 0(x, t < ti) has well defined angular momentum I. Then 0(x, t > ti) 
has well defined angular pseudomomentum m determined by the equation: 



I — m = kn., 

with k E Ij, \m\ < | for even n and \m\ < for odd n. 

Some examples of transformations satisfying this rule have been presented in 
Refs. [32|1M] . Let us analyze an example of a transformation from the point 
of view of angular momentum and pseudomomentum during the evolution. 

The function at t = ti is 0(x, t < ti) = cr^e ''o e*'^ with c = 0.424 and 
To = 2.1. Therefore it has angular momentum I = 3. The function V"(x, t > ti) 



is a squared lattice as in section 3.1, So, since m = I — kn, k G Z, n = A, 



and \m\ < |, necessarily m = —1. In Fig. [3] we present the numerically 
calculated value of / and m during propagation, using Eq. ^ normalized by 
/ dxdyipcj)* to calculate / and discrete rotational transformations to calculate 
m. As can be seen, / for t > ti is no longer conserved but m is. If the definition 
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Fig. 3. Evolution of the angular momentum and pseudomomentum when / = 3 and 
Q is Civ from ti = 1.5. 

of jz is applied to a function with well defined angular pseudomomentum 
0(r, 9, t) = e'^"^^Um{r, 9, t), with Um{r, 9 + e,t) = Um{r, 9, t) the following result 
is obtained: 



So the conservation of m appears as a balance between and jj*. Note that 
J™ is a function of time. 



solution belongs to the corresponding irreducible representation and has the 
proper angular pseudomomentum m = — 1. 



5 Theory of Angular Pseudomomentum in the Quantum Limit 

Once the relevance of a proper treatment of the angular pseudo-momentum 
in the description of discrete symmetry systems has been clearly established, 
one question remains still open: are these properties still valid in the limit 
where the mean-field approach is no longer appropriate? In other words, what 
would happen in the energy domain where quantum effects start to be rele- 
vant? Such a question, that at first glance could seem to be a purely academic 
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one, has become extremely important since Bose-Einstein Condensates (BEC) 
became experimentally feasible thirteen years ago . In order to provide an 
answer to this question one should start by specifying a Hamiltonian to de- 
scribe the interactions at the bosonic level. To do so we have considered a 
standard interaction, the Bose-Hubbard Hamiltonian [38]. This model is a 
discretization of the full quantum many-body Hamiltonian for weakly inter- 
acting ultracold bosons that in the mean-field limit allow us to connect with 
the Gross-Pitaevskii (nonlinear Schrodinger) equation [32] ■ A complete anal- 
ysis of the solution of this model within a periodic ID bosonic chain has 
been carried out [IQ]. From these results it can be clearly observed how the 
aforementioned angular pseudo-momentum relations also hold in the case of 
few-boson systems once the angular pseudo-momentum is properly identified 
with a quantum number by means of the generator of the symmetry group. 
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6 Conclusions 

In conclusion, we have proved a number of theorems and propositions that 
show that the symmetric solitons of a generalized nonlinear Schrodinger equa- 
tion can be classified according to the irreducible representations of the dis- 
crete rotational group of the modulus of their amplitude. We have define 
a new quantity, the angular pseudomomentum, that can be used to classify 
nonlinear solutions in different representations. In this way, the behaviour of 
any nonlinear symmetric sohition under the discrete rotations of the group 
is completely determined by this number. Moreover, we have demonstrated 
several theorems about the allowed values of the angular pseudomomentum. 
We have illustrated these mathematical results with numerical examples in 
and Cqv We have also extended these results to non-stationary solutions 
and we have proved a theorem that states that the angular pseudomomentum 
is dynamically conserved. We have demonstrated how angular momentum and 
pseudomomentum can be related in a (9(2) — C„ interface. Wc have explained 
that the non-conservation of angular momentum in Cnv groups appears as a 
balance between the angular pseudomomentum and another time-dependent 
quantity related to the angularly periodic part {u{r, 6, z)) of the solution. 
Wc have presented a numerical example to illustrate these dynamical results. 
Finally, we have outlined the generalization of these results to the quantum 
limit which is relevant to trapped atoms in an optical lattice when a very low 
number of atoms is considered. It is expected that the most important results 
of this paper survive in this limit since they are preserved by symmetry ar- 
guments that should be impervious to dynamical details, but this remains an 
open question at present. 

LDC acknowledges support by the National Science Foundation under Grant 
PHY-0547845 as part of the NSF CAREER program. 
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